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TECHNICAL FIELD 

This invention relates to graphics, and more particularly to improved 
triangle management in triangular meshes based on a tri-edge structure. 

BACKGROUND OF THE INVENTION 

Computer technology is continually advancing, resuhing in a continuing 
stream of computers that are more powerful than their predecessors. One such 
area of advancement, and one that is of great importance to many designers as 
well as end users, is that of computer graphics. Computer graphics are used in a 
wide variety of fields, including entertainment (e.g., games), computer aided 
design, system modeling, and so forth. 

Typically, computers manage graphics by manipulating small uniform 
geographic shapes that often times are triangles. Graphical objects and surfaces 
are described using groups of these triangles, and the triangles themselves are 
typically organized as a collection of vertices (i.e., each triangle is represented by 
a set of three vertices). Such a representation for triangles is well suited to 
hardware-accelerated rendering. However, such a representation also makes 
determining triangle and vertex adjacency cumbersome and computationally 
expensive. For example, given a particular vertex it can be very time-consuming 
to determine what its edge sharing neighbor vertices are. Determining triangle and 
vertex adjacency is crucial to many algorithms that deal with surfaces (e.g., 
refinement and smoothing algorithms, simplification and level of detail 
algorithms, etc.), yet the current systems for representing triangles are not well- 
suited to determining such triangle and vertex adjacency. 
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The invention described below addresses these disadvantages, providing 
improved structures and processes to manage triangles. 

SUMMARY OF THE INVENTION 

Improved triangle management in triangular meshes based on a tri-edge 
structure is described herein. 

According to one aspect, a data structure having two fields is used to store 
data for each triangle in a triangular mesh. The first field is a set of three vertices 
for the triangle and the second field is a set of three edges, each edge 
corresponding to one of the three vertices. Each of the three edges is an 
identification of a next or subsequent edge that is encountered when performing a 
traversal (e.g., in a counterclockwise direction) about the corresponding vertex. In 
one implementation, each vertex in the set of three vertices includes a set of values 
representing the location of the vertex and an identification of a representative 
triangle edge corresponding to the vertex. 

According to another aspect, triangles can be added to the triangular mesh. 
A triangle to be added to the triangular mesh is identified by a set of three vertices 
for the triangle. Upon receiving the three vertices, a check is made as to whether 
any of the edges that will be part of the new triangle already exist in the triangular 
mesh. If any of these edges already exist, a check is made as to whether 
connectivity of the mesh needs to be changed to accommodate the new triangle, 
and such changes are made if necessary. After performing these checks and any 
necessary connectivity changes are made, additional edges for the triangle are 
created as needed until two edges of the triangle exist. A third edge is then added 
to connect the free ends of the existing two edges, completing the triangle. 
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According to another aspect, triangles can be removed from the triangular 
mesh. The triangle to be removed is indicated by a triangle identifier. For each 
vertex of this triangle, a representative edge for the vertex is updated so that the 
representative edge for the vertex is not an edge of the triangle being removed. 
Then, for each edge of the triangle, the edge is removed and the connectivity of 
the other triangles in the triangular mesh is changed as needed so that any other 
triangle in the mesh having an identification of a next edge that is an edge of the 
triangle being removed has that identification changed to another edge of the 
mesh. 

According to another aspect, three operators are defined to assist in 
management of the triangular mesh: a make edge operator, a splice operator, and 
a swap operator. The make edge operator receives two vertices as inputs and 
generates two triangles, with the two triangles having adjacent edges between the 
two vertices. The splice operator receives two edges as inputs and alters 
connectivity of triangles in the mesh including the edges. The swap operator 
receives a particular edge of a triangle in the mesh as an input, and retums an 
opposite diagonal of a quadrilateral corresponding to the input edge. The make 
edge, splice, and swap operators are selectively invoked to both add triangles to 
the triangular mesh and remove triangles fi-om the triangular mesh. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention is illustrated by way of example and not limitation in 
the figures of the accompanying drawings. The same numbers are used 
throughout the figures to reference like components and/or features. 

Fig. 1 illustrates an exemplary graphics management environment. 
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Fig. 2 illustrates an exemplary triangle and its corresponding tri-edges. 
Fig. 3 illustrates an exemplary structure for storing the data representing a 
triangular mesh. 

Fig. 4 illustrates an exemplary set of operators for traversing edges and 
accessing elements of a mesh. 

Fig. 5 illustrates an example of two triangles having been created. 

Fig. 6 illustrates an example of the operation of a splice operator. 

Fig. 7 illustrates an example of the operation of a swap operator. 

Fig. 8 illustrates an exemplary process for adding a triangle to a triangular 

mesh. 

Figs. 9a and 9b illustrate an exemplary process for checking whether edges 
of a triangle exist. 

Fig. 10 illustrates an exemplary process for determining and changing 
connectivity if needed and recording selected edges. 

Fig. 11 illustrates an exemplary situation that can be encountered with 
respect to the process of Fig. 10. 

Figs. 12a and 12b illustrate an exemplary process for removing a triangle 
from a triangular mesh. 

Fig. 13 illustrates an example of splicing in a new edge during the process 
of Figs. 12a and 12b. 

Fig. 14 illustrates example results of removing edges of a triangle from a 

mesh. 

Fig. 15 illustrates an exemplary computer environment, which can be used 
to implement the processes described herein. 
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DETAILED DESCRIPTION 

Fig. 1 illustrates an exemplary graphics management environment 100 
including an application 102, one or more triangular meshes 104, and mesh 
management modules 106. In one embodiment, meshes 104 are each a set of one 
or more data structures defining the mesh, and are stored in volatile and/or non- 
volatile memory. Application 102 represents any of a wide variety of 
conventional applications, such as games, computer-aided design programs, 
drawing and other graphics programs, etc. Application 102 generates and 
manipulates meshes 104 with the assistance of mesh management modules 106. 

Mesh management modules 106 include an AddTriangle module 108 for 
adding triangles to a mesh 104, a RemoveTriangle module 110 for removing 
triangles from a mesh 104, and one or more operator modules 112 that provide 
various mappings and perform various lower level operations on the triangles in a 
mesh 104. 

Each of meshes 104 is a triangular mesh, which is a set or group of 
triangles that collectively form one or more surfaces when displayed. The 
surfaces can be of any shape - the "triangular mesh" refers to a mesh made up of 
triangles rather than the surface described by the mesh having a triangular shape. 
Each triangle is made up of both a set of vertices and a set of tri-edge structures, as 
discussed in more detail below. Mesh management modules 106 assist in the 
managing of adjacency, also described in more detail below. 

Reference is made herein to homeomorphism. A homeomorphism is a 
continuous 1-1 mapping between spaces whose inverse is also continuous. A 
surface, or 2-manifold, is a point set such that every point has a neighborhood that 
is homeomorphic to a plane. 
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In a triangular mesh, each triangle ^ is a pair of triples, as follows: 

The first triple is a set of three vertices, referred to as a 2-simplex. Each one of the 
vertices v is a pair as follows: 

The value {x,;^} is a tuple of scalar values, such as point coordinates, a normal 
vector, and so forth, and r is a representative edge (this can be any edge that 
begins at the vertex v). 

The second triple in a triangle ^ is a tri-edge structure including three edges, 
each edge being a pair as follows: 

e = {tj} 

where t is the triangle that the edge belongs to, and i is an index value indicating 
the position of the edge in the triangle t. The position index / is an element of the 
set {0,1,2} and mathematical operations described herein involving the position 
index / are taken modulo 3. The edges of a tri-edge structure are not the edges of 
the associated 2-simplex. Rather, each edge is the next edge that would be 
encountered in a counter-clockwise traversal from the associated 2-simplex 
(beginning at the edge of the 2-simplex connecting the vertices v/ and V/_;) about 
the vertex v/. 

In an altemative embodiment, an additional "flip" field / is added to each 
edge, making it a triple rather than a pair (e.g., {t,i/})^ The flip field is a binary 
field to enable the representation of non-orientable manifolds. If the value in the 
flip field is changed, the edges can be viewed as being in the opposite direction 
(e.g., as if the mesh were being viewed from the opposite side). 
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Fig. 2 illustrates an exemplary triangle and its corresponding tri-edges. A 
triangle tk is illustrated as including an edge e = {tj} . As used herein, edges point 
from a vertex V/ (the origin) to a vertex Vt.j (the destination). The tri-edges for 
triangle tj, are also illustrated, with the three edges in the tri-edge structure being 
identified as ei, e^.i, and The edge et (an edge of triangle /;,) is the next edge 
that would be encountered in a counter-clockwise traversal from edge e about the 
vertex v^. Similarly, the edge e^,; (an edge of triangle ty) is the next edge that 
would be encountered in a counter-clockwise traversal from the edge of triangle tj, 
connecting vertices v/_/ and Vi+j about the vertex v/.;, and the edge (an edge of 
triangle Q is the next edge that would be encountered in a counter-clockwise 
traversal from the edge of triangle 4 connecting vertices V/+; and V/ about the 
vertex V/+7. 

A special boundary vertex (also referred to as v^;,/) is used for triangles 
that are on the boundary of the mesh. The boundary of the mesh refers to areas 
where the surface being described by the mesh ends, and can be in the interior or 
at the periphery of the surface. For example, triangles adjacent a hole cut in the 
surface are boundary triangles. Triangles that are not where the mesh ends are 
referred to as interior triangles. Some surfaces may not have any boundary 
triangles (e.g., a sphere), although since such surfaces are typically created 
triangle-by-triangle, such surfaces will have boundary triangles during their 
creation. Any triangle that is at the boundary of a mesh has a vertex of Voo - 
triangles that are on the interior of the mesh do not have any vertex of Voo. In one 
implementation, the value of Voo is a reserved value that a vertex typically does not 
have (e.g., the hexadecimal value ffffffff). 
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Fig. 3 illustrates an exemplary structure 120 for storing the data 
representing a triangular mesh. The data is stored in three portions or lists: a 
vertex list 122, a simplex list 124, and a tri-edge list 126. Vertex list 122 is a list 
of vertices v in the mesh. Simplex list 124 is a list of 2-simplexes in the mesh, 
with each entry in simplex list 124 identifying one of three vertices of the 2- 
simplex. In the illustrated example, each entry in simplex list 124 includes three 
pointers to three vertices in vertex list 122. 

Tri-edge list 126 is a list of tri-edges. Each entry in list 126 is a tri-edge for 
a particular triangle, including three references to {t,i} pairs. The reference to the 
triangle {t) is a pointer to one of the 2-simplexes in list 124, and the position index 
(0 is the numerical value for the position index. A triangle 128 in the mesh 
includes two references - a reference to one of the 2-simplexes in simplex list 124, 
and a reference to one of the tri-edges in tri-edge list 126. 

The structure 120 is merely an exemplary structure for maintaining the data 
representing the triangular meshes. Altematively, any of a wide variety of other 
structures may be used. For example, simplex list 124 may include the actual 
values for the vertices rather than pointers to entries in vertex list 122, By way of 
another example, other structures besides the vertex list, simplex list, and tri-edge 
list may be used. 

Various operators (e.g., operator modules 112 of Fig. 1) are defined and 
used for traversing edges and accessing elements of a mesh. Fig. 4 illustrates an 
exemplary set of such operators. These operators are described with reference to 
an edge e of the triangle having a position index A first set of these operators 
match edges to edges, and are referred to as: rot(e), rot"\e), onext(e), oprev(e), 
and sym(e). 



Iee@hayes piic 509-324-9256 



8 



MS J -731 US PAT APP 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 



The rotate or Yot(e) operator maps to the next edge in the triangle tjc in the 
counterclockwise direction, and is defined as follows: 

The inverse rotate or rof^(e) operator maps to the next edge in the triangle 
tx in the clockwise direction (which is equivalent to twice rotating to the next edge 
in the counterclockwise direction), and is defined as follows: 

vot'\e) = rot(rot(e)) 

The next edge or onext(^) operator maps to the next edge (also referred to 
as the subsequent edge) from edge e when rotating about the vertex V/ in the 
counterclockwise direction. The onext(e) operator is defined as follows: 

onext(e) = tri-edge ei of triangle t^^ 

The previous edge or oprev(£) operator maps to the next edge in the next 
triangle when rotating about the vertex v/ in the clockwise direction. The oprev(e) 
operator is defined as follows: 

oprev(e) ^ rot(onext(rot(e))) 

The same or sym(^) operator maps to the edge of another triangle that 
shares the same vertices as edge e (although the edge is in the opposite direction of 
edge e), and is defined as follows: 

sym(e) = rot(onext(e)) 

Another set of these operators maps edges to vertices. These edge to vertex 
mapping operators are referred to as: org(e), dest(^), right(e), and left(e). 

The origin or org(e) operator maps to the vertex that is the origin of the 
edge e, and is defined as follows: 

org(e) = vertex / of triangle t^^ 
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The destination or dest(6) operator maps to the vertex that is the destination 
of the edge e, and is defined as follows: 
dest(e) = org(rof 

The right(^) operator maps to the other vertex of the edge of triangle t^, that 
shares vertex i with edge e, and is defined as follows: 
right(e) = org(rot(e)) 

The left(e) operator maps to the vertex that is the destination of the next 
edge from edge e when rotating about the vertex V/ in the counter-clockwise 
direction, and is defined as follows: 

left(e) = dest(onext(e)) 

Another set of these operators perform other useful operations, and are 
referred to as: getRep(v), setRep(v,e), setOrg(e,v), and setOnext(a,6). 

The get representative edge or getRep(v) operator maps to the 
representative edge (r) of the vertex v. The set representative edge or setRep(v,e) 
operator sets the representative edge (r) of the vertex v to the edge e. The set 
origin or setOrg(e,v) operator sets the vertex / of the edge e (connecting vertex v/ 
to vertex v,_;) to v. The set next edge or setOnext(a,Z?) operator sets the tri-edge a 
to b. 

Three additional operators are also defined to manipulate triangles in the 
mesh. These three operators are: MakeEdge(vo.V7), Splice(<^,Z?), and Swap(e). 

The MakeEdge(V(?,v/) operator constructs a mesh Tof two triangles a and b, 
and retums an edge of the triangle a. An exemplary implementation of the 
MakeEdge(vo,v;) operator is shown in the following pseudo-code: 
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MakeEdge(vftV7) 

{ 

b ^ {{vo,v;,v«,},{{fl,l},{«,0},{a,2}}}; 
insert a and b into mesh T; 
return {<3,1}; 

} 

The two triangles created by the MakeEdge(v^,v;) operator are illustrated in 
Fig. 5. The MakeEdgeCv^vy) operator creates two triangles a and b having 
adjacent edges between vertex 140 (v^) and vertex 142 (v;). Additionally, each of 
the two triangles a and b share a vertex 144 (Vo^). Although vertex 144 (v^o) is a 
common vertex, it has been shown twice due to the limitations of illustrating the 
triangles on the two-dimensional nature surface of the drawing page. The vertex 
Voo is the special boundary vertex, so the triangle edges are illustrated with dashed 
lines. 

The Splice(a,Z?) operator receives as input a pair of edges a and b and 
rearranges the edge links within the associated triangles so as to re-identify the 
edge pairs. The Splice(a,Z?) operator alters the connectivity of the mesh, changing 
the tri-edge structures for the triangles involved. An exemplary implementation of 
the Splice(a,Z?) operator is shown in the following pseudo-code: 

Splice(a,Z?) 

{ 

a' <— onext(a); 

onext(&); 
a ^ rot(a'); 

/?^rot(Z>0; 
a'^ rot"^(a); 
;5'.-rot-'(^); 

setOnext(a,Z?0; 
setOnext(Z?,a'); 
setOnext(a,y5'); 
setOnext05,a); 
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The effect of the Splice(a,Z?) operator on a pair of simple meshes is 
illustrated in Fig. 6. Given the pair of edges a and b, the Splice(a,Z?) operator 
alters the connectivity of the mesh as illustrated in Fig. 6. The original 
connectivity of the mesh is illustrated by the solid curved lines, while the dashed 
lines indicate new links replacing old links with the same origin. Thus, link 160 is 
replaced by link 162, link 164 is replaced by link 166, link 168 is replaced by link 
170, and link 172 is replaced by link 174. 

In certain embodiments, care should be taken when using the Splice(a,Z>) 
operator to ensure that the cycle of triangles about org(a) and dest(Z?) (that is, the 
triangles encountered when traversing about the vertices org(a) and dest(6)) do not 
overlap, and that the cycle of triangles about dest(a) and org(Z?) do not overlap. 
This can be accomplished at a higher programming level (e.g., by application 102 
of Fig. 1, which calls the Splice(a,6) operator), or alternatively additional checks 
could be added to the beginning of the Splice(a,&) process to verify that these 
cycles of triangles do not overlap (and the process not be performed if they do 
overlap). 

The Swap(e) operator sets the input edge e to the opposite diagonal of a 
quadrilateral and returns this newly set edge. An exemplary implementation of the 
Swap(e) operator is shown in the following pseudo-code: 

Swap(e) 

{ 

a <— onext(^); 

b ^ rot'\e); 

setOrg(a,dest(Z?)); 

setOrg(Z?,dest(a)); 

SpUce(a,^); 

Sphce(^,sym(Z?)); 

retum a; 
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} 

The effect of the Swap(e) operator is illustrated in Fig. 7. The edge 190 (e) 
is the input to the Swap(e) operator. The edge e is then changed and set to be edge 
192 by the Swap(e) operator. 

Given the operators and structures defined above, triangles can be added to 
and removed from a mesh using an AddTriangle process and a RemoveTriangle 
process, respectively. These processes will be discussed in the following figures. 

Fig. 8 illustrates an exemplary process for adding a triangle to a triangular 
mesh. The process of Fig. 8 can be implemented in hardware, software, firmware, 
or combinations thereof 

Initially, a set of three vertices for the triangle to be added are received (act 
220). A check is then made as to whether any of the three edges for the triangle to 
be added already exist (act 222). One or more of the three edges may already 
exist, depending on the vertices for the new triangle and any previous triangles 
added to the mesh. Based on the edges found, a determination is then made as to 
whether the connectivity of any of the pre-existing edges needs to be changed (act 
224). Situations can arise where the triangle to be added cannot be added unless 
the connectivity of some pre-existing edges are changed, as discussed in more 
detail below. In such situations, the necessary changes in connectivity are made 
(act 226). 

A set of zero or more selected edges is then recorded (act 228). These are 
edges that will be used to create the new triangle, as discussed in more detail 
below. The triangle addition process obtains two edges for the new triangle, then 
adds the third edge to those two edges to form the triangle. If two such edges 
already exist, then no new edges need to be created. However, if no edges exist, 
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then two edges need to be created, and if only one edge exists, then one new edge 
needs to be created (act 230). These two edges of the triangle will form a wedge 
or "V" shape. The free ends of the two edges (i.e., the vertex of each edge that is 
not shared by the other edge) are then connected with a third edge (act 232). 

A check is then made as to whether connecting the two edges in act 232 
resulted in filling a triangle-shaped hole (act 234). If the triangle is being added to 
fill in a hole in the mesh of the same size as the triangle being added, then 
additional acts are performed. These additional acts comprise destroying any 
additional edges that were created during the connecting process of act 234 (act 
236). For example, use of the MakeEdge operator discussed above may result in 
having an additional edge (with one vertex that is the boundary vertex) that is not 
needed. After destroying any such additional edges, or if connecting the two 
edges did not result in filling a triangle-shaped hole, then a representative edge is 
set for each vertex of the triangle (act 238), completing the triangle addition 
process. 

Figs. 9a and 9b illustrate an exemplary process for checking whether edges 
of a triangle exist (act 222 of Fig. 8). The process of Figs. 9a and 9b can be 
implemented in hardware, software, firmware, or combinations thereof. The 
process of Figs. 9a and 9b is repeated for each of the three vertices received in act 
220 of Fig. 8 (i.e., the value / ranges from 0 to 2). The process of Figs. 9a and 9b, 
as well as the process of Fig. 10 below, refers to edges a, 6, and c. The edges a, b, 
and c are referred to in Figs. 9a, 9b, and 10 are defined as follows: edge a is the 
edge from vertex to vertex v^+j, edge b is the edge from vertex V/ to vertex V/.;, 
and edge c is the edge from vertex Vf to vertex Voo. 
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Initially, the stored representative edge for the vertex is retrieved (act 
250), and a check made as to whether the representative edge is the empty set (act 
252), If the representative edge is empty (that is, the vertex is an isolated vertex 
that is not part of any edge in the mesh), then a value gi is set to be equal to the 
empty set (act 254), and the process ends for that vertex (act 256). 

However, if the representative edge is not empty, then the representative 
edge is used as a selected edge (act 258) and a variable w is set to be equal to the 
destination vertex of the selected edge (act 260). A check is then made as to 
whether there is an edge from w to the vertex v/_; (act 262). If there is such an 
edge, then a check is made as to whether the vertex to the right in the triangle 
(e.g., using the right() operator discussed above) is equal to Voo (act 264). If the 
vertex to the right in the triangle is not equal to Voo then an error is reported (act 
266) and the triangle addition process ends. However, if the vertex to the right in 
the triangle is equal to Voo then the selected edge is set to be edge b for the new 
triangle (act 268) and a check is made as to whether there are any additional edges 
to select (act 270 of Fig. 9b), In the illustrated example, selection of additional 
edges continues until either all edges have been selected or both of the edges a and 
b have been identified. If there are no additional edges to select then the process 
ends (act 272); otherwise, another edge is selected (act 274) and processing 
continues at act 260 with the newly selected edge. 

Returning to act 262 in Fig. 9a, if there is not an edge from w to the vertex 
V/.; then a check is made as to whether there is an edge from w to the vertex Vi+j 
(act 276). If there is such an edge, then a check is made as to whether the vertex 
to the left in the triangle (e.g., using the left() operator discussed above) is equal to 
Voo (act 278). If the vertex to the left in the triangle is not equal to Voo then an error 
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is reported (act 266) and the triangle addition process ends. However, if the vertex 
to the left in the triangle is equal to Voo then the selected edge is set to be edge a for 
the new triangle (act 280) and a check is made as to whether there are any 
additional edges to select (act 270 of Fig. 9b). 

Returning to act 276 of Fig. 9a, if there is not an edge from w to the vertex 
then a check is made as to whether the destination of the selected edge is equal 
to Voo (act 282 of Fig. 9b). If the destination of the selected edge is equal to v^o then 
the selected edge is set to be edge c for the new triangle (act 284) and a check is 
made as to whether there are any additional edges to select (act 270). Retuming to 
act 282, if the destination of the selected edge is no equal to Voo then a check is 
made as to whether there are any additional edges to select (act 270). 

Fig. 10 illustrates an exemplary process for determining and changing 
connectivity if needed and recording selected edges (acts 224, 226, and 228 of 
Fig. 8). The process of Fig. 10 can be implemented in hardware, software, 
firmware, or combinations thereof The process of Fig. 10 is repeated for each of 
the three vertices received in act 220 of Fig. 8. 

Initially, a check is made as to whether a and b are both not empty (act 
300). If a and b are both not empty, then a check is made as to whether another 
triangle(s) exists between a and b (act 302). Fig. 11 illustrates an exemplary 
situation where another triangle(s) exists between the two edges a and b. As 
illustrated, the edge 310 (a) and edge 312 (Z?) are two edges of the triangle to be 
added to a mesh based on the vertices 314 (v/), 316 (v^.y), and 318 (v/+;). The 
triangles existing in a counter-clockwise traversal about vertex 314 (v^) between 
line 320 (onext(a)) and line 322 (oprev(Z?)) are removed prior to completing the 
triangle addition process. 
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Returning to Fig. 10, if one or more other triangles do exist between a and 
b (act 334) then connectivity of the triangles is changed so that the triangle(s) no 
longer exist between a and b (act 336). A value of gi is then set to be the edge b 
(act 338) and the process ends for vertex V/ (act 340). Returning to act 334, if one 
or more other triangles do not exist between a and b, then processing proceeds to 
act 338 without altering the connectivity of any of the triangles. 

Returning to act 300, if both a and b are not both not empty (that is, one or 
both of a and b is empty), then a check is made as to whether just a is not empty 
and b is empty (act 342). If a is not empty and b is empty, then a value of gi is set 
to be the edge onext(a) (act 344) and the process ends for vertex (act 340). 
However, if a is empty or b is not empty, then a check is made as to whether b is 
not empty(act 346). If b is not empty, then a value of gf is set to be the edge b (act 
348) and the process ends for vertex v/ (act 340). However, if b is empty, then a 
check is made as to whether edge c is not empty(act 350). If c is not empty, then a 
value of gi is set to be the edge c (act 352) and the process ends for vertex V/ (act 
340). However, if c is empty, then an error is reported (act 354) and the triangle 
addition process ends. 

The following pseudocode is an exemplary implementation for the process 
of adding a triangle to a mesh. The pseudocode relies on three predicates, which 
are defined as follows: 

bound(e) = e is not equal to 0 

interior {e) = dest(e) is not equal to Voo 

boundary (e) = dest(e) is equal to Voo 
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AddTriangle(v^),v;,V2) 
{ 

for(/^0,..., 2){ 

e ^ getRep(vO; 
if (bound(e)) { 

do{ 

TV dest(e); 

if (w = V/_;) { 

if (right(^) 7^ Voo) 

Error(" invalid edge"); 

} 

else if (w = v,+;) { 

if(left(e)#v») 

Error("invalid edge"); 

a e; 

} 

else if (w = Voo) e\ 
e ^ onext(e); 
} while {e f eg and not (bound(a) and bound(Z>))); 

gi ^ 0; 

if (bound(a)) { 

if(bound(Z?)) { 

if (onext(^) ^ oprev(Z?)) { 
e ^ onext(^); 

while {{e ^ a) and (dest(e) 7^ Voo)) 
e ^ onext(e); 

if (e = a) 

Error("non-manifold vertex"); 
oprev(Z?); 

Splice(onext(a)j/); 

Splice(e^; 

} 

gi ^ b\ 

} 

else gi ^ onext(a); 

} 

else if (bound(6)) gi ^ b; 
else if (bound(c)) gi ^ c; 
else ErTor("non-manifold vertex"); 
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} 



repeat { 

if (not(interior(g/)) and not(interior(g,./))) { 
e MakeEdge(v„v,./); 
if (boundary(g,)) Splice(onext(e),g/); 
if (boundary^-.;)) Splice(rot"\e),g,.;); 
gi-i ^ TOt-\e); 
gi ^ e; 

} 

else if (interior(g/) and interior(g/./)) { 

Swap(rot"Vg/)); 
if(interiorfe+7)) { 

Splice(rof^fe),g,+/); 

DestroyEdge(g,+;); 

} 

setRep(v,,g/); 
setRep(v,+;,rot(g,)); 
setRep(v,+2,ror'(g,)); 
return; 

} 

i <—/+!; 



Figs, 12a and 12b illustrate an exemplary process for removing a triangle 
from a triangular mesh. The process of Figs. 12a and 12b can be implemented in 
hardware, software, firmware, or combinations thereof 

Initially, an identifier of the triangle to be removed is received (act 402). 
For each vertex in the identified triangle, the representative edge for the vertex is 
updated to be either an edge that is not part of the identified triangle or the empty 
set (act 404), so that once removed no vertex will have a representative edge that 
is an edge of the removed triangle (and thus no longer part of the mesh). Also, for 
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each vertex V/ in the identified triangle a temporary variable {hf in the illustrated 
example) is set up holding the value of oprev(e/) for the edge (act 406). 

A check is then made as to whether the identified triangle is completely 
surrounded by other triangles (act 408). If so, then a new edge is spliced into the 
mesh having a third vertex of Voo (act 410), allowing new triangles to be created 
(having a vertex at Voo) as part of the removal process. Fig. 13 illustrates the 
process in act 410 in more detail. Assume that a triangle is going to be removed 
from the interior of mesh 420 (a cross-hatch pattern is used in Fig. 13 to identify 
the triangle to be removed). This will leave a hole in mesh 420 where the triangle 
was, and so an additional edge is added as illustrated in mesh 422 having a third 
vertex at Voo. This creates the basis for three additional triangles in the hole being 
left by the triangle being removed, and each of the three additional triangles is a 
boundary triangle. 

Returning to Fig. 12a, after splicing in the new edge in act 410, or if the 
identified triangle is not completely surrounded by other triangles, then one vertex 
of the identified triangle is selected (act 430). A check is then made as to whether 
right(/z/) given the selected vertex equals Voo (act 432). If right(/z/) given the 
selected vertex does not equal Voo then another vertex is selected (act 434), and the 
check is repeated. Once a vertex is selected such that right(/?/) equals Voo then the 
Swap operator (discussed above) is used to remove the edge hi from the mesh (act 
436). This situation is illustrated in more detail in Fig. 14. The value of right(/2/) 
where hi is edge 430 is the vertex 432. By using the Swap operator, the edge 430 
is removed and replaced with the edge 434 (having a vertex at Voo). 

Returning to Fig. 12b, after removing hi a check is made as to whether 
right(/2y-,;) equals Voo (act 438). If right(A/+;) equals Voo then triangles with the edge 
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hi^i are removed (act 440) and the edge hi+j is destroyed (act 442). Referring to 
Fig. 14, the resultant mesh after removing the two triangles that had edge hi+i as 
an edge and destroying edge hi-^i is mesh 446. 

Retuming again to Fig. 12b, after destroying edge in act 442 (or if 
right(/?/+7) does not equal Voo), then a check is made as to whether right(A/.7) equals 
Voo (act 448). If right(/2/_/) does not equal Voo then the removal process ends (act 
450). However, if right(/^/_;) equals Voo then triangles with the edge A/.; are 
removed (act 452), the edge is destroyed (act 454), and the removal process 
ends (act 450). Referring to Fig. 14, the resultant mesh after removing the two 
triangles that had edge hi^ as an edge and destroying edge is mesh 456 (just the 
three vertices v^, V/_;, and Vy+7. 

The following pseudocode is an exemplary implementation for the process 
of removing a triangle fi*om a mesh. 

RemoveTriangle(0 
{ 

for(z = 0,...,2){ 

hi <— oprev(^); 
V,- ^ org(^); 
if(right(/zO = Voo) 

if(oprev(/z/) = onext(e))) 
setRep(v/,0); 

else 

setRep(v/,oprev(oprev(/z/))); 

else 

setRep(V;,A,); 

} 

if ((right(Ao) ^ Vco) and 
(right(/z;)^ Voo)and 
(right(/2,)7^vj) { 

e ^ MakeEdge(v(),V7); 

Splice(/zo,e); 
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} 

repeat { 

if(right(AO^v«,){ 
Swap(/z/); 

if (right(/z/+7) = Voo) { 
e ^ sym(//K;); 
Splice(oprev(/z/+;),rot"^(e)); 
Splice(oprev(e),onext(^)); 
DestroyEdge(/z/+;); 

} 

if (right(/z/.7) = Voo) { 
e ^ sym(hi.i); 
Splice(oprev(A,-.7),rof^(e)); 
Splice(oprev(e),onext(e)); 
DestroyEdge(/z/.7); 

} 

return; 

} 

} 



The manifold spaces represented using the structures, operators, and 
processes described herein exhibit the following characteristics: 

a. The triangle mesh resulting from a series of calls to MakeEdge() and 
valid calls to Splice() is homeomorphic to a collection of closed, oriented surfaces. 

b. The triangle mesh resulting from a series of calls to AddTriangle() 
and RemoveTriangleO is homeomorphic to a collection of closed, oriented 
surfaces such that all vertices are shared by at least three triangles (although an 
exception to this is that AddTriangle() will allow the creation of a simple mesh 
component with all three edges shared by a pair of interior triangles). 
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c. The triangle mesh resulting from a series of calls to AddTriangle() 
and RemoveTriangleO with boundary triangles removed is homeomorphic to a 
collection of oriented surfaces, possibly with boundary. 

The structures described herein can be implemented in any of a wide 
variety of manners. The following are exemplary implementations of the 
structures using the C++ programming language, although other structures could 
altematively be used. The following is the structure definition for a mesh 
structure: 

struct Mesh 

{ 



Vertex[] 


vertices; 


Edge[] 


representatives; 


Triangle[] 


triangles; 


Triedge[] 


triedges 


uint32 


numVertices; 


uint32 


numTriangles; 


uint32 


numlnterior; 


uint32 


numBoundary; 



{; 

In the Mesh structure definition, "uint32" refers to an unsigned 32-bit 
integer, "Edge[]" is a set of edge structures (defined below) for the mesh, 
'Triangle[]" is a set of triangle structures (defined below) for the mesh, and 
'Triedge[]" is a set of tri-edge structures (defined below) for the mesh. "Vertex[]" 
is a set of vertex structures for the mesh, each identifying a vertex in the mesh (the 
definition of the vertex structure is system dependent, and thus is not discussed 
further herein). Additionally, "numVertices" refers to the total number of vertices 
in the mesh, "numlnterior" refers to the number of interior triangles in the mesh, 
"numBoundary" refers to the number of boundary triangles in the mesh, and 
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"numTriangles" refers to the total number of triangles in the mesh (should be equal 
to the sum of numlnterior and numBoundary), 

The "Edge[]'\ "Triangle[]'\ and "Triedge[]" structures are defined as 
follows: 

struct Edge 
{ 

unsigned int t : 30; 
unsigned int i : 2; 

{; 



Struct Triangle 

{ 

uint32 vert[3]; 

{; 



struct Triedge 
{ 

Edge[] edge[3]; 



Thus, as can be seen by the definitions, each edge includes a triangle 
identifier t and a position index each triangle contains three vertex structures, 
and each tri-edge contains three edge structures. 

Fig. 15 illustrates an exemplary computer environment 500, which can be 
used to implement the processes described herein. The computer environment 500 
is only one example of a computing environment and is not intended to suggest 
any limitation as to the scope of use or functionality of the computer and network 
architectures. Neither should the computer environment 500 be interpreted as 
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having any dependency or requirement relating to any one or combination of 
components illustrated in the exemplary computer environment 500. 

Computer environment 500 includes a general-purpose computing device in 
the form of a computer 502. Computer 502 can be used to implement, for 
example, environment 100 of Fig. 1. The components of computer 502 can 
include, but are not limited to, one or more processors or processing units 504, a 
system memory 506, and a system bus 508 that couples various system 
components including the processor 504 to the system memory 506. 

The system bus 508 represents one or more of any of several types of bus 
structures, including a memory bus or memory controller, a peripheral bus, an 
accelerated graphics port, and a processor or local bus using any of a variety of 
bus architectures. By way of example, such architectures can include an Industry 
Standard Architecture (ISA) bus, a Micro Channel Architecture (MCA) bus, an 
Enhanced ISA (EISA) bus, a Video Electronics Standards Association (VESA) 
local bus, and a Peripheral Component Interconnects (PCI) bus also known as a 
Mezzanine bus. 

Computer 502 typically includes a variety of computer readable media. 
Such media can be any available media that is accessible by computer 502 and 
includes both volatile and non-volatile media, removable and non-removable 
media. 

The system memory 506 includes computer readable media in the form of 
volatile memory, such as random access memory (RAM) 510, and/or non- volatile 
memory, such as read only memory (ROM) 512. A basic input/output system 
(BIOS) 514, containing the basic routines that help to transfer information 
between elements within computer 502, such as during start-up, is stored in ROM 
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512. RAM 510 typically contains data and/or program modules that are 
immediately accessible to and/or presently operated on by the processing unit 504. 

Computer 502 may also include other removable/non-removable, 
volatile/non- volatile computer storage media. By way of example, Fig. 15 
illustrates a hard disk drive 516 for reading from and writing to a non-removable, 
non-volatile magnetic media (not shown), a magnetic disk drive 518 for reading 
from and writing to a removable, non-volatile magnetic disk 520 (e.g., a "floppy 
disk"), and an optical disk drive 522 for reading from and/or writing to a 
removable, non- volatile optical disk 524 such as a CD-ROM, DVD-ROM, or other 
optical media. The hard disk drive 516, magnetic disk drive 518, and optical disk 
drive 522 are each connected to the system bus 508 by one or more data media 
interfaces 526. Alternatively, the hard disk drive 516, magnetic disk drive 518, 
and optical disk drive 522 can be connected to the system bus 508 by one or more 
interfaces (not shown). 

The disk drives and their associated computer-readable media provide non- 
volatile storage of computer readable instructions, data structures, program 
modules, and other data for computer 502. Although the example illustrates a hard 
disk 516, a removable magnetic disk 520, and a removable optical disk 524, it is to 
be appreciated that other types of computer readable media which can store data 
that is accessible by a computer, such as magnetic cassettes or other magnetic 
storage devices, flash memory cards, CD-ROM, digital versatile disks (DVD) or 
other optical storage, random access memories (RAM), read only memories 
(ROM), electrically erasable programmable read-only memory (EEPROM), and 
the like, can also be utilized to implement the exemplary computing system and 
environment. 
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Any number of program modules can be stored on the hard disk 516, 
magnetic disk 520, optical disk 524, ROM 512, and/or RAM 510, including by 
way of example, an operating system 526, one or more application programs 528, 
other program modules 530, and program data 532. Each of such operating 
system 526, one or more application programs 528, other program modules 530, 
and program data 532 (or some combination thereof) may implement all or part of 
the resident components that support the distributed file system. 

A user can enter commands and information into computer 502 via input 
devices such as a keyboard 534 and a pointing device 536 (e.g., a "mouse"). 
Other input devices 538 (not shown specifically) may include a microphone, 
joystick, game pad, satellite dish, serial port, scanner, and/or the like. These and 
other input devices are connected to the processing unit 504 via input/output 
interfaces 540 that are coupled to the system bus 508, but may be connected by 
other interface and bus structures, such as a parallel port, game port, or a universal 
serial bus (USB). 

A monitor or other type of display device 542 can also be connected to the 
system bus 508 via an interface, such as a video adapter 544. In addition to the 
monitor 542, other output peripheral devices can include components such as 
speakers (not shown) and a printer 546 which can be connected to computer 502 
via the input/output interfaces 540. 

Computer 502 can operate in a networked environment using logical 
connections to one or more remote computers, such as a remote computing device 
548. By way of example, the remote computing device 548 can be a personal 
computer, portable computer, a server, a router, a network computer, a peer device 
or other common network node, and the like. The remote computing device 548 is 
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illustrated as a portable computer that can include many or all of the elements and 
features described herein relative to computer 502. 

Logical connections between computer 502 and the remote computer 548 
are depicted as a local area network (LAN) 550 and a general wide area network 
(WAN) 552. Such networking environments are commonplace in offices, 
enterprise-wide computer networks, intranets, and the Internet. 

When implemented in a LAN networking environment, the computer 502 is 
connected to a local network 550 via a network interface or adapter 554. When 
implemented in a WAN networking environment, the computer 502 typically 
includes a modem 556 or other means for establishing communications over the 
wide network 552. The modem 556, which can be intemal or external to computer 
502, can be connected to the system bus 508 via the input/output interfaces 540 or 
other appropriate mechanisms. It is to be appreciated that the illustrated network 
connections are exemplary and that other means of establishing communication 
link(s) between the computers 502 and 548 can be employed. 

In a networked environment, such as that illustrated with computing 
environment 500, program modules depicted relative to the computer 502, or 
portions thereof, may be stored in a remote memory storage device. By way of 
example, remote application programs 558 reside on a memory device of remote 
computer 548. For purposes of illustration, application programs and other 
executable program components such as the operating system are illustrated herein 
as discrete blocks, although it is recognized that such programs and components 
reside at various times in different storage components of the computing device 
502, and are executed by the data processor(s) of the computer. 
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Computer 502 typically includes at least some fomi of computer readable 
media. Computer readable media can be any available media that can be accessed 
by computer 502. By way of example, and not limitation, computer readable 
media may comprise computer storage media and communication media. 
Computer storage media includes volatile and nonvolatile, removable and non- 
removable media implemented in any method or technology for storage of 
information such as computer readable instructions, data structures, program 
modules or other data. Computer storage media includes, but is not limited to, 
RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, 
digital versatile disks (DVD) or other optical storage, magnetic cassettes, magnetic 
tape, magnetic disk storage or other magnetic storage devices, or any other media 
which can be used to store the desired information and which can be accessed by 
computer 502. Communication media typically embodies computer readable 
instructions, data structures, program modules or other data in a modulated data 
signal such as a carrier wave or other transport mechanism and includes any 
information delivery media. The term "modulated data signal" means a signal that 
has one or more of its characteristics set or changed in such a manner as to encode 
information in the signal. By way of example, and not limitation, communication 
media includes wired media such as wired network or direct-wired connection, 
and wireless media such as acoustic, RF, infrared and other wireless media. 
Combinations of any of the above should also be included within the scope of 
computer readable media. 

The invention has been described herein in part in the general context of 
computer-executable instructions, such as program modules, executed by one or 
more computers or other devices. Generally, program modules include routines, 
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programs, objects, components, data structures, etc. that perform particular tasks 
or implement particular abstract data types. Typically the functionality of the 
program modules may be combined or distributed as desired in various 
embodiments. 

For purposes of illustration, programs and other executable program 
components such as the operating system are illustrated herein as discrete blocks, 
although it is recognized that such programs and components reside at various 
times in different storage components of the computer, and are executed by the 
data processor(s) of the computer. 

Ahematively, the invention may be implemented in hardware or a 
combination of hardware, software, and/or firmware. For example, one or more 
application specific integrated circuits (ASICs) could be designed or programmed 
to carry out the invention. 

The discussions herein describe various operators and exemplary 
implementations that are specific to certain directions. In alternate embodiments, 
these directions can be different. For example, the Onext() operator is discussed 
with reference to counterclockwise traversal about a vertex. Altematively, the 
OnextQ operator may refer to clockwise traversal about a vertex. 

Conclusion 

Although the description above uses language that is specific to stmctural 
features and/or methodological acts, it is to be understood that the invention 
defined in the appended claims is not limited to the specific features or acts 
described. Rather, the specific features and acts are disclosed as exemplary forms 
of implementing the invention. 
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